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This paper studies the inter-annual precipitation variations in different regions of East Asia from oceans to 
interior areas in China during 1979 - 2012. The results computed by Empirical Orthogonal Functions (EOF) 
demonstrate that the annual precipitation changes are mainly related to the El Nino-Southern Oscillation, 
East Asian summer monsoon and aerosols. We also found that the increased Sea surface temperature (SST) 
could explain the precipitation changes over the Northwest Pacific in the dry season (Oct. - May) and the 
East China Sea and the South China Sea in the rainy season (Jun. - Sep.). The precipitation changes over the 
ocean unexplained by SST were likely due to the water vapor transport dominated by dynamic factors. With 
the increased SST, the moisture transported from oceans to interior land was likely redistributed and caused 
the complicated regional variability of precipitation. Moreover, the impacts of aerosols on cloud and 
precipitation varied with different pollution levels and different seasons. 

Ocean precipitation accounts for approximately 80% of global precipitation and plays a vital role in the 
Earth's climate system 1 . Its variations can directly impact the intensity of the global hydrological cycle 
and the redistribution of precipitation over land and sea 2 . Enormous efforts have been made to improve 
the understanding of how ocean precipitation responds to warming climate 36 . Globally, precipitation increases in 
equatorial rain bands; decreases in subtropics as greater tropical convection in the rising branch of the Hadley 
circulation will lead to enhanced subsidence in the subtropics; and increases in high-latitudes due to increase in 
moisture transport 7-5 . However, the response of precipitation to the global warming varies regionally. For the 
tropical ocean, two viewpoints have coexisted in interpreting precipitation changes, i.e., one predicts intensified 
rainfall in presently rainy regions (wet-get-wetter) 3,6 , and the other suggests that increased precipitation would 
occur where the sea surface temperature (SST) rises more (warmer-get- wetter) 4,5 . Huang et al. 7 recently proposed 
that the two kinds of mechanisms above are not contradictory, but complementary, i.e., the variability of the 
annual mean precipitation over the tropical ocean follows warmer-get-wetter, whereas the seasonal-mean pre- 
cipitation follows the wet-get-wetter mechanism. Moreover, the precipitation in the mid- and high latitudes has 
increased by about 0.5 to 1% per decade, except over East Asia 8 . What causes the exception in East Asia is still 
poorly understood 10,11 . 

East Asia is a region experiencing huge human population and rapid economic growth in recent years. The 
precipitation in East Asia can be affected by several climatic factors. El Nino-Southern Oscillation (ENSO), a 
coupled ocean-atmosphere tropical Pacific phenomenon, is known to affect regional precipitation by altering the 
tropical cyclones in different ways 12 . As the strongest signal of inter-annual climate variability, ENSO may be 
becoming more frequent in the warming climate (cycle from 7-8 years to 3-4 years) 13 . In addition, the monsoons, 
i.e., East Asian summer monsoon (EASM), the South Asian summer monsoon (SASM), the Indian summer 
monsoon (ISM) and the western North Pacific summer monsoon (WNPSM), are also cyclical factors greatly 
affecting regional precipitation in Asia and all have a seasonal cycle. Recent report indicates that the EASM has 
experienced significant weakening after 1970s, thus leading to a decrease of precipitation in north China 14 . The 
SASM is also weakening and is found to be influenced by the increase of anthropogenic aerosols 15 . Unlike ENSO, 
EASM and SASM etc., SST experienced a monotonic increase. The increase is proposed to have a strong control 
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on the spatial variation of precipitation 13 , and the relationship 
between the warming SST and precipitation is nonlinear 16 . 

The anomalous changes of regional precipitation under global 
warming may not be solely due to the climate variabilities such as 
ENSO, monsoon and SST changes. Atmospheric aerosols can also 
play a role 15 , which are highly varying with time and space. They can 
impact regional precipitation by affecting the cloud microphysics 
properties and even weakening the monsoons 15 . Moreover, atmo- 
spheric aerosols may partially offset the positive forcing due to green- 
house gases 17 . Asia, particularly China, has become one of the highest 
aerosol loading regions in the world. The increased summer floods in 
south China and droughts in north China have been associated with 
the increasing loadings of aerosols 18 . 

The regional precipitation has a direct influence on agriculture 
activities and local economic and social development. Large uncer- 
tainties are still present on the forecast of regional precipitation even 
for the most advanced models 11,20 , although the responses of precip- 
itation to ocean warming has been extensively studied on the global 
scales 5,19 . An in-depth analysis of long-term variability of regional 
precipitation in East Asia and related climatic and environmental 
factors are crucial to reveal the intrinsic mechanisms and would 
benefit for forecasting the future. In view of this, the objectives of 
this study are to characterize inter-annual variability of precipitation 
over various monsoon-influenced zones in East Asia during 1979- 
2012 and to improve the understanding of inhomogeneity of regional 
precipitation changes related to ocean warming. In addition, we will 
pay particular attention to the influences of the cyclical factors 
(ENSO and the monsoon circulations), aperiodic climate factor 
(SST) and regional factors (aerosols and clouds) on the changes of 
precipitation over our study zone in the last three decades. 

Methodology 

Data sources. The datasets used in this study include: 1) version 2.2 
monthly average Global Precipitation Climatology Project (GPCP) 
products (available at http://precip.gsfc.nasa.gov/) on a 2.5° X 2.5° 
grid. GPCP is the combination of a variety of micro-wave data, 
infrared data and the surface rain-gauge data 21 . GPCP preci- 
pitation data has been extensively used for the precipitation 
analysis in the tropics and mid-latitudes 21,22 . 2) ERSST (version 3b) 
which is the US National Oceanic and Atmospheric Administra- 
tion's (NOAA) Extended Reconstruction of SST dataset (http:// 
www.esrl.noaa.gov/psd/), with a 2.0° X 2.0° resolution 23 . The 
ERSST analysis is based on the available Global Telecommuni- 
cations System (GTS) ship and buoy observations from the Inter- 
national Comprehensive Ocean- Atmosphere Data Set (ICO ADS). 
3) NOAA Nino-3.4 index which is defined as the average SST 
anomalies over the Nifio-3.4 region (5°S-5°N, 170°E-120°W). The 
SST variability in the Nino-3.4 region is reported to have a strong 
effect on shifting rainfall in the western Pacific and is often used as an 
ENSO index 24 . 4) Four monsoon indices are used in this study to 
examine the relationship between the precipitation and monsoons, 
such as EASMI, SASMI, IMI and WNPMI. The definition of the four 
monsoon indices has been listed in Table 1. 5) Monthly average 
Moderate-resolution Imaging Spectroradiometer (MODIS) Level 3, 
Collection 5 datasets onboard the Terra satellite (http://ladsweb. 
nascom.nasa.gov/data/search.html). It has a spatial resolution of 



1.0° X 1.0°. The specific parameters used in this study include 
aerosol optical depth (AOD) at 550 nm, fine mode (aerodynamic 
diameter (Dp) < 1 urn) AOD and clouds fraction (CF), cloud 
effective radius (CER). MODIS aerosol products can monitor the 
temporal and spatial characteristics of the global aerosol field over 
land and ocean. In particular, the MODIS aerosol products over 
China have been extensively validated using the Aerosol Robotic 
Network (AERONET) 29 and the Chinese Sun Haze meter Network 
(CSHNET) data 30 . 6) Global reanalysis data sets provided by 
National Centers for Environmental Prediction/National Center 
for Atmospheric Research(NCEP/NCAR) for the period 1979- 
2012 are used to estimate the terms of the atmospheric water 
vapor transport 31 . The physical variables used in this study include 
the monthly average specific humidity, and the meridional and zonal 
wind components at eight standard pressure levels: 1000, 925, 850, 
700, 600, 500, 400, and 300 hPa. The vertically-integrated moisture 
flux can be expressed as 



Q=- 



q Vdp, 



and the moisture flux divergence can be represented as 
V(Vq/g)=^(Vq/g), 



(1) 



(2) 



where Q is the vertically integrated moisture flux, q is specific 
humidity, V is the horizontal wind vector, p is the pressure, p s is 
the surface pressure, g is the acceleration due to gravity, and 
Vp*( Vq/ g) is the divergence of the integrated moisture flux. 

Study areas. Our study regions include both the oceans and the 
continents in East Asia (Fig. 1). The oceanic zones contain the 
Yellow Sea (YS, 120°E-126°E, 32°N-39°N), East China Sea (ECS, 
117°E-129°E, 25°N-34°N), South China Sea (SCS, 106°E-121°E, 
3°N-22°N) and Northwest Pacific (NP, 122°E-147°E, 1°N-35°N). 
Over continents, we focus on the coastal areas within 250 km from 
YS, ECS and SCS. Over the coastal areas the summer moisture is 
mainly from the ocean, especially the coastal areas of ECS and SCS, 
but the loadings of atmospheric aerosols are much higher compared 
to those in the marine atmosphere. For inter-comparison, the 
precipitation changes over the monsoon-influenced regions in 
China are also analyzed. These regions are interior lands, but the 
moisture is partially influenced by the monsoons. The interior lands 
are further divided into north China and south China according to 
the traditional climate regionalization. The Qinling-Huaihe Line is 
the demarcation line between south and north China 32 . 

Methods 

Empirical orthogonal functions (EOF) and the Pearson product-moment correlation 
analysis (PCC) are used for data interpretation. The EOF is a statistical method that 
can extract the dominant modes of change from a time-dependent two-dimensional 
dataset. The EOF analysis gives spatial patterns of each mode and the time series of the 
loadings of the modes. The EOF spatial patterns are eigenvectors, which are ortho- 
gonal so that the leading EOFs describe the spatially coherent pattern that maximizes 
its variance 33 . The PCC is usually used to measure the degree of linear relationship 
between two variables 34 . The correlation will be noted when the statistical significance 
is equal to or exceeds the 95% confidence level. 



Table 1 | The definition of four monsoon indices used in this study 
Monsoon Indices 



Definition 



Time domain 



Reference 



East Asian Summer Monsoon Index (EASMI) U850(l 1 0°-140°E,l 0°-40°N) JunetoAugust 25,26 

South Asian Summer Monsoon Index (SASMI) U850(35 C -97.5°E,5 -22.5 C N) June to September 25,26 

Indian Monsoon Index (IMI) U850(40 C -80°E,5 3 -1 5"N)-U850(70°-90 o E,20 o -30°N) June to September 27,28 

Western North Pacific Monsoon Index (WNPMI) U850(l 00°- 1 30°E,5°-1 5°N)-U850(1 1 0 -1 40°E,20°-30 :, N) June to September 27,28 

Note: U850 is defined as the area-averaged dynamical normalized seasonality (DNS) zonal wind anomalies at 850 hPa. 
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Figure 1 | Map of the study regions (the grid points for precipitation data reading are marked in different colors, representing different areas) . The map 

in the figure was created by D.Y. using Arcgis 9.3. 



To measure the atmospheric water-holding capacity as a function of local T, we use 
an approximation to the Clausius-Clapeyron relation of the form: 

e s =Aexp[BT/(C + T)], (3) 

where e s is the saturated vapor pressure (kPa), T is the surface air temperature( L 'C) 
and the constants A, B, and C have the values 0.6108, 17.27, and 237.3, respectively 35 . 
Since marine air temperature has a good agreement with SST in daily average over the 
majority of large ocean regions, SST data can be used as a proxy for marine air 
temperature. This is especially true for this study as the time resolution of SST used in 
this study is monthly mean 36 . We compute the increasing rate of the calculated water- 
holding capacity per °C SST (AW/ AT) using the formula: 

AW 

^yjT — ( e s2012 ~~ e sl979)/ e sl979/C^2012 — T 1979 ), (4) 

e sl979 and e s2012 indicate the e s in 1979 and 2012, T 1979 and T 2012 stand for the SST in 
1979 and 2012, respectively. The net increase of precipitation by the net increase of 
SST (AP/AT) is computed with the formula: 

AP 

^ ~(P2012 — ^1979)/^1979/C^2012 — ^1979)' w) 

Pi 979 and P 2012 are the precipitation in 1979 and 2012. The units of AW/ AT and AP/ 
AT are both %/°C. Globally, this equation determines that the water-holding capacity 
of the atmosphere increases by about 7% for every 1°C rise in temperature from 1900 
to 2002 37 . 

Results and Discussion 

Regional variability of precipitation in East Asia related to cyclical 
factors. To explore temporal and spatial distributions of annual 
precipitation variations in East Asia during 1979-2012, EOF 
analysis is performed on annual precipitation anomaly fields. Only 
the first three EOFs and corresponding principal components (PCs) 
(Fig. 2) have passed the North test 33 and they account for 52.8% of the 



variance. The three EOFs and corresponding PCs are discussed in 
terms of correlations between the PCs of each mode and the possible 
influencing factors (Nino-3.4 index, EASMI, SASMI, ISMI, 
WNPSMI and AOD, CER, CF). The first EOF accounts for 33% of 
the anomaly variance (Fig. 2a) and shows a negative phase over the 
NP, the SCS, the north YS and parts of the north China. The variation 
of the precipitation is the product of the space field and time field of 
the EOF. The corresponding PCI (Fig. 2d) has a decreasing tendency 
from 1979 to 2012, implying that the precipitation in these regions 
has increased in the past three decades. The increase is more obvious 
over the NP and SCS than over the north YS and parts of the north 
China. A significant correlation (r = 0.53) is obtained between the 
PCI and the Nino-3.4 index from non-trend interannual variations, 
indicating that the first mode largely reflects ENSO-related trend 
variations. Smith et al. 38 and Trenberth et al. 39 have analyzed the 
global precipitation using EOF and also attributed the first mode 
to the influence of ENSO on the basis of correlations between the 
PCs and the Nino index. 

The second mode (Fig. 2b) seems to be associated with the ENSO 
as well 38,39 , except that PC2 shows a significant negative correlation (r 
= —0.58) with the Nino-3.4 index. Furthermore, PC2 also shows a 
negative correlation to the WNPSMI (r = -0.62) and the EASMI (r 
= —0.4), but has a significant positive correlation with CF (r = 0.64). 
There is a weaker East Asian monsoon circulation and a delayed 
onset with the warmer SST anomalies in the region of Nino-3.4 40 . 
Thus, the indices are not independent from each other and all have a 
significant correlation with PC2. Excluding the two ENSO-related 
abnormal years, i.e., 1999 (flood year) and 2005 (drought year), the 
PC2 shows a slight decrease before 1995, but then a slight increase. 



SCIENTIFIC REPORTS | 4:5693 | DOI: 1 0. 1 038/srep05693 



3 



first EOF of East Asia (33S) 





SO 1985 11)90 1995 2000 2005 2010 



The second EOF of East Asia (13. 6N) 






80 1985 1990 1995 2000 2(M)5 2010 



The third EOF of East Asia (6 2\) 

"if"" ' ^ 



/ 






f 




0}» 




8 


«> 




6 


on 










1 


ut 










2 


0<S 






1 




> 0 


-oos 




-2 






-1 


-0»> 










-6 


■05 






( M 




ft 




I9S5 1 990 1995 2000 2005 2010 



Figure 2 | The first three EOFs (left panel, (a), (b), (c)) and principal components (right panel, (d), (e), (f)) for annual precipitation of study areas 
during 1979-2012, and the red line in the figure indicates the linear fitting of the principal components. The maps in the figure were created by D.Y. 
using Matlab 2011b. 



The decreasing trend in the southern part of the NP in the second 
mode could be partly associated with the reduced water vapor 
amount owing to the weakening of EASM, especially the 
WNPSM 14 . There is an eastward-movement of the westernmost 
588-dagpm contours of the 500 hPa geo-potential height field over 
western Pacific (approximately 134°E before 1993 and 132°E after 
1994) (also see Supplementary Fig. SI), which indicates that the 
WNPSM is weakening. Although the water vapor amount also 



reduced over the ECS, SCS and their coastal areas (Figs. 3 and 4), 
the precipitation over these areas got enhanced. That can occur 
because the precipitation is determined by both local evaporation 
amounts and the moisture fluxes from remote moisture source 
regions. 

The third mode (Fig. 2c) accounts for 6.2% of the total variance 
and the corresponding PC3 (Fig. 2f) exhibits an obvious increasing 
trend. No significant correlations are found between PC3 and Nino- 
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Figure 3 | The vertically integrated mean water vapor transport over our study areas in the dry season (Oct.-May). (a), the vertically integrated 
mean water vapor transport for 1979-1993. (b), the difference of the vertically integrated mean water vapor transport between 1979-1993 and 1998-2012. 
The arrow key is standing for the water vapor flux and the color shading is representing the water vapor flux divergence. The unit is g/(cnrhPa - s) for the 
water vapor flux and g/(cm 2- hPa - s) for the water vapor flux divergence. The turning time chosen for the comparison is the shift time of the SST over the 
four seas, i.e., 1993-1998 in the dry season. The maps in the figure were created by G.T. using Matlab 201 lb. 

3.4 index or between PC3 and EASMI, but PC3 has a significant are also reported in EOF analysis of global precipitation by Gu and 
correlation with total AOD (r = —0.64). This suggests that the third Adler et al. 41 . In which, two modes appears to be related to ENSO, 
mode is possibly related to the atmospheric aerosols. Three modes similar to our study. However, they attributed the third EOFs to some 
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Figure 4 | Same as Fig. 3 but for the rainy season (June-Sep.), and the turning time is 1994. The maps in the figure were created by G.T. using Matlab 
2011b. 



aspects of regional precipitation changes. The third EOFs in this 
study shows a stronger increase in NP, but a weaker decrease in 
the south China and YS. We will discuss the relationships between 
precipitation and atmospheric aerosols in a later section. 

Regional variability of ocean precipitation related to SST. For 

oceans, the elevated SST is the most direct indicator of global 



warming. Climate models and satellite observations both indicate 
that the total amount of water in the atmosphere increased at a 
rate of 7% per kelvin of surface warming and the observed 
precipitation also exhibited an increasing rate of 6% per kelvin of 
surface warming over the past two decades 19 . Westra et al. 42 reported 
that the median intensity of extreme precipitation changed 
significantly in proportion to changes in global mean temperature 
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at a rate between 5.9% and 7.7% per kelvin. The larger moisture 
holding capacity of air in a warmer climate can cause more 
precipitation. With increasing global air temperature, the 
thermodynamic component associated with an increase of water 
vapor generally shows a positive contribution to precipitation 
changes 6 . However, the dynamic component due to the change in 
atmospheric motion could reduce or increase precipitation, 
depending on locations 7-9-43 . 

We first examine the variability of SST and then examine possible 
relationship between SST and precipitation in the dry and rainy 
seasons separately. The homogeneity tests can find out the qualitative 
change of the dataset in the time series analysis by computing the 
variances 44 . An obvious shift during 1993-1998 detected by the 
homogeneity tests occurs both in the dry and rainy seasons (Fig. 5) 
indicates that the substantive change of the precipitation over the 
four oceans take place in 1993-1998. This shift appears to be assoc- 
iated with a climate regime shift in the Pacific basin. Bond et al. 45 
reported a similar SST shift in 1998/99 over the western and central 



NP which was characterized by a meridional seesaw of SST anom- 
alies. In the dry season, the shift over the YS, ECS, SCS and NP occurs 
in 1993/1994, 1996/1997, 1996/1997 and 1997/1998, respectively. 
For rainy season, the shifts all occur in 1993/1994 regardless of the 
locations of our study oceans. The shifts occur earlier in oceanic 
regions are closer to the continent. The precise reason for this shift 
remains unclear. After the SST jump, a downward trend is observed 
for the dry season SST at a rate of 0.22-0. 33°C per decade over our 
oceanic regions, whereas in the rainy season SST fluctuates with 
no clear trend. The shift of the SST in the dry or rainy season may 
have an influence on the precipitation, even though they are not 
coincident. 

To further study the impact of the SST on precipitation, we divided 
the net increase of precipitation by the net increase of SST (AP/AT) 
and compared with the increasing rate of the calculated water-hold- 
ing capacity per °C SST (AW/ AT). In the dry season, the ratios of AP/ 
AT are —0.02 ± 0.003%/°C (mean ± standard deviation) over the 
YS, - 14 ± 2%/°C over the ECS, 75 ± 4.4%/°C over the SCS and 23 ± 




26 0 
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Figure 5 | Comparison between the dry season (Oct.-May) precipitation and SST over Yellow Sea, East China Sea, South China Sea, the 250 km 
coastal areas adjoining to three seas, the south China, the north China, and Northwest Pacific during 1979-2012. The red dot line stands for the linear 
trend of the precipitation. The red line in the figure indicates the mean of the SST or precipitation before the shift tested by the homogeneity test with the 
Excel 2010, while the blue line indicates the means of the SST or precipitation after the shift. 
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1%/°C over the NP. The mean and standard deviation were calcu- 
lated from the ratios over seven sub regions of the YS, 13 sub regions 
of ECS, 34 sub regions over the SCS and 88 sub regions over the NP 
(Fig. 1). The rates of AW/ AT are 10 ± 1.4%/°C over the YS, 15 ± 
6.9%/°C over the ECS, 22% ± 4.3/°C over the SCS and 21 ± 0.7%/°C 
over the NP, respectively. For the NP, the mean value of AP/AT is 
close to the mean value of AW/ AT, indicating that the precipitation 
changes over the NP in the dry season is dominated by the SST. Adler 
et al. 46 also show a positive correlation between the SST and precip- 
itation over tropical oceans. However, the mean value of AP/AT over 
the SCS is approximately 3.5 times of the mean value of AW/AT. 
Dynamic components may enhance the effect of warming SST on 
precipitation over the SCS (Table 3). Moreover, more precipitation 
over the SCS in the dry season may be related to the Subtropical high 
of the NP moving westward (see Supplementary Fig. S2). The intens- 
ive moisture in the dry season over the SCS is actually originated 
from the subtropical ocean, see also in Fig. 3. On the contrary, the 
mean value of AP/AT over the YS and ECS are negative. Dynamic 
components, e.g., circulation changes in a small regional scale, appar- 
ently overwhelm the local thermodynamic components to reduce the 
precipitation over the two regions. In the rainy season, over the ECS, 
the mean value of AP/AT (19 ± 0.9%/°C) is almost the same as the 
mean value of AW/ AT (21 ± 0.8%/°C). The same is true for SCS. We 
postulated that precipitation changes over the ECS and SCS in the 
rainy season are probably determined mainly by thermodynamic 
components. The mean value of AP/AT is 42 ± 1.4%/°C over the 
NP, which almost doubles the corresponding mean value of AW/AT. 
The increase in inflowing moisture from the southeast along the 
southwest flank of the West Pacific Subtropical high likely caused 
the additional increase of AP/AT ratio relative to AW/ AT ratio (see 
Fig. 4). Over the YS, the mean value of AP/AT is -5.7 ± 0.9%/°C, 
countrary to the mean value of AW/ AT (17 ± 2.8%/°C). Again, the 
weakening mid-latitude westerly may dominate the moisture trans- 
port over the YS (Fig. 4) and lead to the exceptions in this region. 

Regional variability of coastal and interior land precipitation 
related to ocean precipitation. The coastal and interior land 
precipitation is directly associated with agriculture activities, social 
and economic development. In China, approximate 30% of 
population resides in coastal areas (within 250 km of the YS, ECS 
and SCS). In the dry season (see in Fig. 5), the precipitation over the 
coastal area adjoining to the YS decreases by 10 ± 0.4 mm per decade 
and the rate is approximately 30% smaller than over the YS. 
However, over the interior zone (>250 km from the YS) and the 
whole monsoon-influenced region of north China, the precipitation 
increases by 9 ± 2.4 mm per decade in the dry season. For the coastal 
area adjoining to the ECS, the precipitation increases by 6 ± 0.8 mm 
per decade in the dry season. The discrepancy between the positive 
rates of AW/ AT and the negative ratios of AP/AT in the dry season 
over the YS and the ECS suggests that less moisture is transported by 
the weakening mid-latitude westerly (see Fig. 3). Moreover, the 
increase of precipitation in north China has also been proposed to 
be due to the Pacific Decadal Oscillation 47 , but a poor correlation 
(0.05) between the precipitation in north China and the Pacific 
Decadal Oscillation obtained in this study does not support this 
hypothesis. For the coastal area adjoining to the SCS and the whole 
monsoon-influenced area in south China, the precipitation between 
has a positive correlation (r = 0.78) in the dry season and decreased 
by 58 ± 4.1 mm per decade over the former and 27 ± 4.3 mm per 
decade over the later area. This suggests that the dry-season moisture 
source for the two regions is primarily the mid-latitude westerly from 
the continent (see Fig. 3) which has weakened in the last three 
decades. Recall that the mean value of AP/AT over the SCS is 
approximately 3.5 times the mean value of AW/ AT in the dry 
season. More precipitation over the SCS might also reduce the 
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moisture from the SCS being transported to inland areas through the 
water vapor passing along the SCS (see Fig. 3). 

In the rainy season (Fig. 6), the coastal area adjoining to the YS has 
a slight increase in precipitation, yet the precipitation over the north 
China decreases by 12 ± 1.3 mm per decade. More droughts during 
summer has been reported in north China since the late 1990s 48 . The 
physical connections, if any, among precipitation decrease over the 
YS, increase over the coastal area adjacent to the YS and decrease over 
the north China need more modeling studies which are beyond the 
scope of this study. The rainy season precipitation of the coastal area 
adjacent to the ECS exhibits an increase with the rate 30 ± 1.2 mm 
per decade and the rate is almost the same as that over the ECS. The 
supply of water vapor from the SCS and the NP to ECS and its coastal 
areas is greatly reduced (Fig. 4). The increase over the ECS is likely 
due to its increasing moisture holding capacity of air as mentioned 
earlier. The same reason could also be true for the increase of pre- 
cipitation over the coastal area adjacent to the ECS. The rainy season 
precipitation over the coastal area adjacent to the SCS has a shift in 
1993/1994 and the shift matches the rainy season SST shift in the 
SCS. The rainy season precipitation over the coastal area from the 
SCS increases by 40 ± 1.7 mm per decade, which is much larger than 



the 12 ± 6.3 mm per decade increase over the SCS in the same 
season. The precipitation over the monsoon-influenced zones in 
south China decreases by — 10 ± 2.5 mm per decade. The moisture 
transported by cross-equatorial flow and South Asian monsoon to 
the SCS were largely weakened in the rainy season (see Fig. 4), also 
supporting that the increase of the precipitation over the SCS may be 
due to the thermodynamic components. Moreover, Nanling 
Mountains can partially block the moisture from the SCS in the rainy 
season, leading to more precipitation in the windward side and less 
precipitation in the leeward side. Overall, with the increasing SST, 
moisture transport from oceans to interior land is likely redistribu- 
ted, causing the complicated variability of precipitation over different 
regions in the dry and rainy seasons. 

Variability of precipitation related to aerosols. Aerosols can act as 
cloud condensation nuclei (CCN) and have an influence on regional 
precipitation 49 . The effects of aerosols on the precipitation is mostly 
to suppress it 9 , but they can sometimes enhance it 49 . Using a 
combination of observations and modelling results, Yi et al. 51 
proposed that there is a negative correlation between AOD and 
precipitation over land and a positive correlation over ocean. The 
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Figure 7 | Time series of the monthly AOD anomaly during 2000-2012 over the study areas in the dry season (Oct.-May), and the red line in the figure 
indicates the linear fitting of the AOD anomaly, and AOD represents the average of the AOD value in different areas. 

AOD data are only available after 2000, and thus our study is they are 0.7 ± 0.2 over the coastal area adjoining the YS, 0.5 ±0.1 
confined to 2000-2012, with results shown in Figs. 7 and 8. In the over the coastal area adjacent to the ECS, 0.5 ± 0.2 over the coastal 
dry season, the values of AOD are 0.5 ± 0.1 over the YS, 0.4 ± 0.1 area adjacent to the SCS, 0.3 ± 0.1 over north China, and 0.4 ± 0.1 
over the ECS, 0.2 ± 0.08 over the SCS, 0.2 ± 0.05 over the NP while over south China, respectively (Figs. 7 and 8). Tang et al. 52 classified 
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Figure 8 | Same as Fig. 7 but for the rainy season (June-Sep.). 
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Table 3 | The correlations between the precipitation and related factors (CF, CER, AOD, Fine-Mode (aerodynamic diameter (Dp) < 1 \im)) 
AOD and SST) from the NP to the interior China in the dry season (Oct.-May) and rainy season (June-Sep.) during 2000-201 2. All the 
correlation coefficients are tested using the Pearson correlation test, the single asterisk is for P < 0.05, and the double asterisk is on behalf of 
P < 0.01 

Dry Season Rainy Season 

Fine-Mode Fine-Mode 
CF CER AOD AOD SST CF CER AOD AOD SST 



PRCP YS 






0.33** 


0.27** 




YS 


0.57** 


0.31* 


0.29* 


ECS 




0.68** 


0.29** 


0.23* 


0.21* 


ECS 


0.73** 




0.31* -0.65* 


SCS 


0.75** 


0.70** 


-0.46** 


-0.33** 


0.51** 


SCS 


0.34* 




0.40** 


NP 


0.48** 


0.75** 


-0.38** 


-0.34** 


0.68** 


NP 


0.72** 






250 km YS 






0.39** 


0.24* 




250 km YS 


0.57** 


0.28* 


0.47** 


250 km ECS 


0.57** 


0.64** 


0.34** 


0.31** 




250 km ECS 


0.74** 


0.30* 


0.38** 


250 km SCS 


0.47** 


0.67** 








250 km SCS 


0.72** 


0.64** 


-0.33* 


north China 


-0.38** 


-0.40** 


0.31** 


0.40** 




north China 


0.79** 




0.41** 


south China 


0.54** 


0.83** 


0.39** 






south China 


0.78** 


0.83** 





AOD into three types: low (0-0.3), moderate (0.3-0.6), and high 
(0.6-0.9) pollution levels. According to the classification, the 
average value of AOD over the SCS, NP, and north China was low 
(Fig. 7). The average value of AOD in the coastal area adjacent to the 
YS was the highest and the area that had the highest emissions of S0 2 
in China 53 . The remaining AODs rank at the moderate level. In the 
rainy season, the values of AOD are 0.5 ± 0.2 over the YS, 0.3 ± 0.08 
over the ECS, 0.2 ± 0.06 over the SCS, 0.8 ± 0.3 over the coastal area 
adjacent to the YS, 0.5 ±0.1 over the coastal area adjacent to the ECS, 
0.5 ±0.1 over the coastal area adjacent to the SCS, 0.3 ±0.1 over 
north China, 0.1 ± 0.04 over the NP, and 0.4 ± 0.07 over south 
China. The ECS, SCS, and NP are all clean areas by this 
classification. When the time series of AOD anomaly is examined, 
AOD anomaly first increased and then decreased over the SCS and 
NP in the dry season during 2000-2012, and the turning point was 
around 2006. No clear temporal trend for AOD anomaly is exhibited 
over other regions. In the rainy season the increasing trend is found 
only over the YS and the coastal area adjacent to the YS. The AOD 
anomaly over the NP increased first and then decreased after 2006, 
similar to the change in the dry season. The other areas exhibit a 
decreasing trend in AOD anomaly. 

To further investigate the role of aerosols on the cloud microphy- 
sical properties and the effects of aerosols and clouds on precipita- 
tion, the correlations among the MODIS satellite AOD, CER and CF 
are computed for 2000-2012 (Table 2) as well as the correlations 
between the precipitation and AOD, CER, CF in the same regions 
(Table 3). From the perspective of aerosols, the atmosphere over the 
SCS and the NP is relatively clean, the coast areas adjacent to the YS 
suffers from severe particulate pollution, while the remaining study 
regions are moderate. In the relatively clean atmosphere, more aero- 
sols are usually accompanied by the smaller CER, less CF and sup- 
pression of precipitation in the dry season. Rosenfeld et al. 49 pointed 
out that the fewer aerosol droplets, combining with the adequate 
moisture caused by the increasing SST, could make CCN grow and 
drop to the ocean surface quickly. The hypothesis can explain the 
negative correlations between AOD and other factors (i.e., precipita- 
tion, CER and CF) over the SCS and NP. For the rainy season, more 
aerosols are generally associated with smaller CER in the clean areas, 
but have no significant correlation with the precipitation. The sum- 
mer monsoons may have a strong influence on the precipitation over 
the SCS and NP in the rainy season, overwhelming the effects of 
aerosols. With the adequate water vapor transported by the summer 
monsoons, greater CF is associated with more precipitation in the 
rainy season (Table 3). For the severe pollution area with high AOD 
loadings, i.e., the coastal area adjoining the YS, more aerosols and 
bigger CER probably promotes precipitation in both the dry and 
rainy season. Traditionally, excessive aerosols combining with the 



increasing water vapor may form a large number of smaller CCN, 
prolonging the lifecycle of the cloud, and reducing the precipita- 
tion 49 . However, Koren et al. 50 recently proposed that excessive aero- 
sols combining with high water vapor can sometimes enhance 
precipitation. For the remaining areas, the relationship among the 
precipitation, AOD, CER and CF are complicated and location- 
dependent. The best correlation of 0.83 between the precipitation 
and CER in the dry season occurs in south China. In the rainy season, 
the AOD has a positive correlation with the precipitation only in ECS 
and the coastal area adjacent to the ECS. The fine mode AOD is also 
used for correlation analysis. In theory, the fine mode aerosols should 
have better correlations with the cloud and precipitation than the 
total aerosols since fine mode particles usually have much higher 
number concentration. However, the correlations between the 
fine-mode AOD and the CER, CF, precipitation are very similar with 
the total AOD in both dry and rainy seasons. Overall, precipitation 
over the continents is likely affected more by aerosol and related CF 
and CER, while precipitation over the oceans is likely affected more 
by climate factors. As shown, the influence of aerosols on the pre- 
cipitation varies with pollution level and season. 

Conclusions 

This work combines statistics analysis of satellite data and diagnostic 
analysis of global reanalysis data sets for water vapor transport to 
study the regional precipitation changes from the Northwest Pacific 
(NP) to the interior of China under a warming climate and discusses 
the possible influence on the precipitation from four factors: ENSO 
(El Nino-Southern Oscillation), monsoons, SST (Sea surface tem- 
perature), and aerosols. Our results show that the first Empirical 
Orthogonal Functions (EOFs) of precipitation is mainly related to 
the ENSO, the second EOFs is associated with the ENSO, East Asian 
Summer Monsoon (EASM), and Western North Pacific Monsoon 
( WNPM), whereas the third EOFs is possibly dominated by the effect 
of aerosols. The role of South Asian Summer Monsoon (SASM) and 
Indian Monsoon (IM) on the precipitation in the monsoon-influ- 
enced zones is largely overwhelmed by the factors mentioned above, 
although the moisture from the Indian Ocean was one of the major 
sources for precipitation in East Asia. 

The increased SST is an important cause for precipitation changes 
over the NP in the dry season and over the East China Sea (ECS) and 
South China Sea (SCS) in the rainy season. Over these oceanic 
regions, the ratios of the net increase of precipitation by the net 
increase of SST (AP/AT) are almost the same as the ratios of the 
increasing rate of the calculated water-holding capacity per °C SST 
(AW/AT). For the SCS in the dry season and the NP in the rainy 
season, the ratios of AP/AT are 2-3 times AW/ AT and dynamic 
components likely contributed more to the increased ratios of API 
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AT. The Subtropical high of the NP moving westward may cause 
additional increment of AP/AT ratios over SCS in the dry season 
while increase in southeast inflowing moisture along the southwest 
flank of the West Pacific Subtropical high may result in the additional 
increment of AP/AT ratios over the NP in the rainy season. However, 
the ratios of AP/AT are negative in both dry and rainy seasons over 
the YS, although the rates of AW/ AT are highest among the four 
oceanic regions, suggesting that dynamic components and other 
factors overwhelm local thermodynamic components in those 
regions. In addition, the rainy season precipitation for coastal areas 
adjacent to the ECS exhibits an increase, with the rate almost same as 
that over the ECS. The rainy season precipitation over the coastal 
area adjacent to SCS has a shift in 1993/1994 and the shift matches 
the rainy season SST shift in the SCS. This evidence suggests that the 
increased SST not only affects ocean precipitation, but also the pre- 
cipitation over the adjacent coastal areas. 

In the relatively clean atmosphere, i.e., over the SCS and the NP, 
larger aerosol optical depth (AOD) are usually associated with the 
smaller cloud effective radius (CER), less clouds fraction (CF) and 
suppression of precipitation in the dry season. In the rainy season, no 
clear association was found for those regions. The coastal area adja- 
cent to the YS is severely pollution area with high AOD loadings. Its 
larger AOD is associated with bigger CER and likely promotes the 
precipitation in both the dry and rainy seasons. For the remaining 
areas, the relationship among the precipitation, AOD, CER and CF 
are less clear. 
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